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Abstract. We study two-dimensional systems with boundary curves described 
by power laws. Using conformal mappings we obtain the correlations at the bulk 
critical point. Three different classes of behaviour are found and explained by scaling 
arguments which also apply to higher dimensions. For an Ising system of parabolic 
shape the behaviour of the order at the tip is also found. 



The shape of a system undergoing a second order phase transition can have a 
strong influence on its critical behaviour. This is shown by the results for edges 
(in 3D) or corners (in 2D). The local critical exponents are then continuous functions 
of the corresponding angle [1-6]. But this striking feature also raises the question 
which property of the boundary actually causes it and what would be obtained for 
other shapes. For the critical behaviour long range effects are essential and thus a 
simple rounding of the corner will not matter [7]. We therefore study here shapes 
which differ from the corner geometry in the large: the boundary curves are described 
by power laws and do not have asymptotes. The prototype is the parabola. We 
use conformal mappings to obtain the critical correlation functions for various two- 
dimensional geometries. When the system forms the interior of a general parabolic 
figure, we find a new unusual form of the critical behaviour. When the system forms 
the exterior, on the other hand, one recovers the behaviour of a system with either a 
straight surface or a cut. These results can be understood from the way the boundary 
curves behave under renormalization. A similar classification will therefore hold in 
three dimensions. Our results at the critical point are complemented by a calculation 
of the tip magnetization for an Ising model of parabolic shape which also shows unusual 
features. 

Consider first a system with free boundaries in the form of a simple parabola 
v 2 = 2pu + p 2 in the plane w = u + iv as in figure (la). It can be related to the 
half-plane z = x + iy,y>0by the conformal map 
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(a) (b) (c) 

Figure 1. The three types of geometries considered in the text. 

At criticality, the correlation function in the half-plane has the form [3] 

G{^z 2 ) = { yi y 2 )- x i)(u) 



(2) 



where the scaling function depends on the variable u> = 4y 1 y 2 /\z 1 — z 2 \ 2 and has 
the asymptotic form ifj((v) ~ u Xl for small u>. Here x and x 1 denote the bulk and 
surface scaling dimensions, respectively, of the operators in G. Using the standard 
transformation [8], one finds in the parabolic geometry 



G(w 1 ,w 2 ) = A^A^ cosh ( x cosh ( 2 cos t] 1 cos rj 2 ) x ifj((v) 
where the rescaling factors are 

7T / sinh Q + sin rj i 



(3) 



A, 



4p V C 2 + r% 



(4) 



and we have used parabolic coordinates n^J ^ = £ + \rj. The system is then 

characterized by < i] < n/2. If w 1 , w 2 lie on the positive w-axis, the variable 10 
is given by u = 4 cosh cosh ( 2 /( cosh ( x — coshC 2 ) 2 . One then finds for ( x fixed and 
C 2 > 1 (u 2 > p) 



G(w 1 ,w 2 ) = A^)—^ exp 



u 



-7TX-, 




(5) 



This is not a simple power law in u 2 as one finds for a corner, which clearly shows the 
difference between the two cases. It also differs from the result in a strip. However, if 
both u l ,u 2 ^> p and ^fu^ — ^Jul ^> v /p, it can be written as 
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with L(u) = 2y/2pu being the width of the system at position u. If in addition 
u 2 — u l <C (u 2 + u 1 )/2 = u, an expansion around u in (6) gives back the strip result [8]. 
In this sense, the parabola can be considered as a strip of varying width. 

If one fixes the boundary variables one can also discuss the order parameter profile. 
A transformation as in (3) then gives, for w = u > 



<*(«)> = A 



7r 2 tanh ( 
4K 



(7) 



Thus ((p) ~ u~ x l 2 for u ^> p and there is no exponential factor in this case. The 
exponent, however, is still different from its value x found for a corner. 

One can easily generalize the treatment to boundaries which have the asymptotic 
form 



±Cu c 



a < 1. 



(8) 



One merely has to shift the parabola into the right half-plane and then to distort it. 
This amount to the replacement 



w 
2p 



w 
2p 



2(l-a) 



1 



(9) 



in equation (1). The quantity C is then C = {2p) 1 a /2(l — a). This changes the 
result (5) for G into 



G{w 1 ,w 2 ) = A( Ul ) — exp 



—nXi 



2C(1 -a) 



(10) 



The functional form of G thus varies continuously with the parameter a describing 
the boundary shape. The more this shape approaches the corner geometry (a — > 1), 
the slower the exponential falloff becomes. On the other hand, for a = the system 
forms a half-strip and one recovers the corresponding simple exponential decay [8,9]. 
For a < 0, it has a spoon-like shape and the decay becomes very rapid in the narrow 
region. 

We now turn to a system in the shape of figure (lb). Here the boundary is curved 
towards the outside, so that a > 1 in equation (8). To relate it to the upper z-plane 
one has to use a different mapping, namely 



z = 1 



w 



1/8 



(11) 



where s = 1 — 1/a. Asymptotically, one now has z = iw and the rescaling factor 
|dw/d^| appearing in the transformation of G becomes one. Therefore, for w 1 ,w 2 on 
the axis, with u 1 fixed and u 2 — > oo, one always obtains the result of the half-plane 



G(w 1 ,w 2 ) = A{u x )- 



U 



(12) 



In this sense, this type of boundary is equivalent to a straight surface. 
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Finally, a system with a cut-out portion in the form of equation (8) as in figure 
(lc), can be obtained via a mapping 



z — 1 



w s 



(l)] 1/2s («) 



where now s = 1 — a and a < 1 again. Asymptotically, the relation therefore is, for 
all a, z = iy/w. But with such a transformation one maps the ^-half-plane onto the 
w-plane with a cut [2-4]. The correlation function therefore is asymptotically 

G( Wl ,w 2 )=A( Ul )-^ (14) 
u 2 

with the corner exponent corresponding to the cut. 

The preceding results can be understood if one considers the behaviour of the 
boundary curve, equation (8), under a change of scale in a renormalization procedure. 
With u' = u/b, v' = v/b it becomes 

v ' = ±b a - 1 C(u') a (15) 

or 

C' = b a ~ x C (16) 

Thus, for a > 1, C grows under renormalization and the boundary curve approaches a 
straight line. For a = 1, C is invariant and thus a marginal variable. This explains the 
particular role of a corner formed by two straight lines. For a < 1, C decreases and 
the system approaches either a cut geometry or a one-dimensional line geometry. In 
the latter case, however, one has a non-ordering system and this causes the particular 
features of the parabolic geometry. 

According to equation (16), 1/C may be considered as a scaling field with 
dimension 1 — a (like 1/L in finite-size scaling). It vanishes at the half-plane fixed 
point. One may therefore write the following scaling ansatz for the correlations along 
the •u-axis 

<»^) = ^ g (t-t-^)- < 17 » 

With b = C 1 /i-« ? on e gets 

«(«,«4)=^ /1 -^(l^).I^)) («) 

where L(u) = 2Cu a is the width of the system at u. Equation (6) can thereby be 
generalized to any value of a < 1 with the scaling function given by 



g(a v a 2 ) ~ ( ai a 2 ) xa/1 a exp 



TTX 1 



(19) 



when a l5 a 2 ^> 1 and a 2 — a l ^> 1. This can also be verified explicitly. The scaling 
behaviour of the order parameter profile is obtained in the same way and reads 

<*(«)> = L{u)-*f (20) 
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where, according to equation (7), lim a ^ 00 f(a) = 0(1). 

Relations like (18) and (20) correspond to a local formulation of finite-size 
scaling. One may also notice that all these scaling considerations still apply in higher 
dimensions. 

Finally, let us address briefly the ordered state. A priori, it is not obvious that 
a system in the shape of figure (la) will order at all. We have therefore studied 
an Ising model with parabolic shape v = ±Cy/u, using the corner transfer matrix 
technique [10]. This means that one considers the transfer matrix connecting the 
spins at the upper and lower boundaries with fixed boundary condition on the right 
end of the system. Assuming a square lattice in the Hamiltonian limit [11] one then 
is lead to study the following operator (describing an inhomogeneous transverse Ising 
chain) 

"JV-l N-l 

h = -c j2^<+ x J2^+ T <<+i ( 21 ) 

. n=l n=0 

where A -1 measures the temperature and N is the size of the system along the axis. 
The coefficients reflect, in the sense of a continuum limit, the number of vertical and 
of horizontal bonds at position n, respectively. The transverse field vanishes at n = 
due to the absence of vertical bonds for the first spin and at n = N as a consequence 
of the boundary condition. The operator can be diagonalized in terms of fermions. 
The single-particle excitation energies t v = 1Cu) v then follow from 

nAV£_i + <^ + 1)C + n\r n+ i = "lV n (22) 

with appropriate boundary conditions at n = 0, N. This system of equations is similar 
to one studied previously in a related context [12] and, as there, can be solved with 
Gottlieb polynomials. In the limit N — > oo one finds, for A > 1 



v = V(A 2 - l)u, z/ = l,2,3--- (23) 



Identifying the boundaries, the magnetization at the tip of a system which is isotropic 
at the critical point, is given by [10,11] 

™ = <O = n tanh (y)- ( 24 ) 

V 

Evaluating this near the critical point (A > 1) leads to 

(25) 



m Q ~ exp 



C 2 (A-1) 



where a = 7£(3)/16 ~ 0.526. Thus there is order, but it vanishes exponentially fast 
at the critical point. This reflects the difficulty to maintain it in such a geometry. We 
note that the argument in the exponential can be expressed as the ratio £/p where 
£ ~ (A — 1) _1 ~ t~ v is the bulk correlation length. 

The behaviour of the tip magnetization may also be deduced from scaling 
considerations. The magnetization at position u along the axis satisfies 
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which, with b = t v ', leads to 



M^W/fe^H (27) 



where (3 and z/ are bulk exponents. 

One may even go further assuming that, at the tip, the leading contribution to 
the magnetization which is induced by the bulk at a distance D ~ (C/C) 1 ^ (where 
the width of the system is of the order of the bulk correlation length) , decays with D 
like the correlation function in (10) when u 1 — ► and u 2 = D ^> C 1 / 1- ". Then 

m (t, 0, 1) = m Q ~ exp (- a Jl—^ (28) 

and the temperature dependence of the tip magnetization follows 

t -u(l-a) \ 1 l° r 



m ~ exp 



—a 



C(l - a) 



(29) 



in agreement with (24) for the Ising parabola with a = | and v = 1. 

One should mention that the results (10) for the correlation function and (29) for 
the order parameter are quite similar to those obtained for an Ising model with bond 
strengths decreasing towards a free surface as K(n) = K(oo)(l — A/n y ), y < 1 [13-15] 
with the correspondences y «-> a, A <-> C. This can be understood qualitatively since 
in both cases the surface order near the critical point can only be maintained through 
the action of the far-away bulk portion of the system. 

Finally, for an anisotropic system with correlation length exponents z/y (along the 
it-axis) 7^ u ± , the scaling dimension of 1/C is changed into 1 — av^/v_ L . Therefore the 
perturbation to the half-plane geometry then becomes relevant when a < v±_jv^. 
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